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ON THE INVERSE PROBLEM OF THE CALCULUS OF 

VARIATIONS * 

By G. A. Bliss 

In the second volume of the Annals, t Prof. Osgood has shown in a very 
interesting way why an arc 

V = y( x )> x £x£x lt 

which gives a minimum (or a maximum) value to an integral of the form 

(1) /?(*» y» y*)dx, 

must be an extremal, that is, a solution of the differential equation 

d 

(2) ?f-as*v = o- 

This equation is of the second order, and its solutions are therefore a two-para- 
meter family of curves. 

In the so-called inverse problem of the calculus of variations the two- 
parameter family of solutions of equation (2) is given in advance, and it is 
required then to find the integrals of the form (1) for which these curves are 
the extremals. In connection with some problems in differential geometry 
DarbouxJ has found that there is an infinity of such integrals, and showed 
how they may be determined. For more general cases the corresponding 
problems are discussed in an interesting article by Hirsch.§ 

In the present paper it is proposed to apply the method of Darboux to 
integrals in the form 

( 3 ) //(*.y. T)\£fTyi<a. 

"The first part of this paper, $§1, 2, was presented before ( The Chicago Section of 
the American Mathematical Society, April 22, 1905, under the title " The inverse problem of 
the calculus of variations In parametric representation." 

tSer. 2, vol. 2, p. XOTj 1901. 

XLecotu sur la thiorie giairale des surface*, vol. 3, §{604-605. 

§ Mathematitehe Annalen, vol. 49 (1897), p. 49. 

(127) 
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The arcs along which the integral is taken are supposed to be represented 
in the form 

(4) * = *(<), y = y(t), toZt&ti, 

and t is the angle between the tangent and the x-axis defined by the 
equations 

«< • Vt 

cos t = , , sin t = 



V^ + yf V^ + yi' 

or, what is the same thing, by 

(5) T = To+ i, 4 + y\ dt > 

where t is the initial value of t at the point t = t . The integral (3) evi- 
dently reduces to (1) when 

/(*» y» t) = ^(*, y, tanr) cost. 

The geometrical advantage in using an integral of the form (3) instead 
of {1) is due to the fact that the arcs over which it is taken may go in any 
direction, while in the form in which Darboux considered the problem diffi- 
culties are always encountered when the tangent to the curve becomes parallel 
to the y-azis, that is, when y x becomes infinite. For the inverse problem the 
integral (3) is also decidedly more convenient than the well-known integral 

(6) f >(x, y, x u y t )dt 

A 

which is usually used when the curves are taken in parametric representation. 
This is because the function f is determined by integrating a single linear 
partial differential equation instead of the system of such equations which 
present themselves when the solution of the problem for the integral (6) is 
attempted. 

The first two sections below are devoted to the derivation for the integral 
(3) of the so-called Euler equation corresponding to equation (2), and to the 
determination of the integrals in the form (3) which have a given family of 
extremals. The integrals which Darboux gives and those which are found in 
§2 involve two arbitrary functions. In §3 it is shown what further conditions, 
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besides the prescription of the extremals, suffice to determine the integral 
uniquely, and in §4 the results of the preceding sections are applied to the 
determination of the integrals which are minimized by straight lines. 

1. Sulei^B equation. The derivation of Euler's equation for the 
extremals of the integral (3) has been made in a previous paper,* but for the 
sake of completeness it will be repeated here. An arc E in the form (4) is 
supposed to give to the integral a smaller valuo than that given by any other 
curve in its neighborhood with the same endpoints. A variation V of E in 
the form 

(V) X=x(t) + a£(t), Y=y(t) + a V (t), t £ t £ t u 

where f and i\ are two functions of t vanishing at t and t lt can be made to lie 
as near as is desired to E by taking the parameter a sufficiently small, and for 
a — the arc V coincides with E itself. It is evident then that the integral 



1(a) = r\f(X, Y, T) \IX\ + Yi dl 



is a function of a which must have a minimum, and consequently the derivative 

-j- must vanish, for a = 0. 
da 

The derivative can be readily calculated when it is noticed that 

-3T« = £ » X u = £<* Y a = t), Y u = ti„ 
and that 

T — J±?? ~ ^'^ 
2 *~ X\ + Yi 

since 

omT^ — S , sinr= T * 



vTf+TT MX) + Y\ 

Then 

[ifo] "/r'* + f" V + (/ COST _ /^ 8inT )f» + (fairiT+f T co8T)t),\d8, 
where 

ds = \Jx> t + yt dt, 

* A generalization of the notion of angle, Transactions of the American Mathematical 
Society, vol. 7 (1906), p. 186. 
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or by a partial integration 

<?> [3...-r{[ / --^ c/co8T - / * 8iDT) ] { 

+ \f' ~ Ts ( ^ 8iD T + /' COS T )] ' } dt ' 

This expression must be zero however f and n are chosen, and by the usual 
argument in the calculus of variations it follows that along the arc E the two 

expressions 

P = / x - -£ (/cos t -/ T sin t) , Q =/„ - ^ (/sin t +/ t cos t), 

must vanish. These conditions upon E are not independent, for the relation 

(8) Pcost + QsinT = 

holds between P and Q, and one verifies readily that the vanishing of the two 
expressions P and Q is equivalent to the vanishing of the expression T de- 
fined by the equation 

(9) Psin-r — Q cost = T. 
In fact from (8) and (9) 

P= Tsinr, ^=-Tcost. 

From either of these equations the value of Tis easily calculated. 

Alowj any arc E which minimizes the integral (3) the expression 

(10) T(x, y, t, t.) =/, sin t -/, cost + f„ cos t + f yr sin t + (f + f„)r, 
must vanish identically. The equation 

(11) T(z, y, t, t.) = 

is called the Elder equation for the integral (3). 

2. Determination of the integrals which have given ex- 
tremals. Consider a family of curves 

(12) x = x(t, u, v), y = y(t, w, v), 
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involving the two parameters u and v. From equation (5) the direction 
angle r is also a function of t, u, v, 

(13) T=\(t, u, v), 

and it will be supposed that the three equations (12) and (13) have single- 
valued solutions 

(14) t = t(x, y, t), u = u(x, y, t), v = v(x, y, t) 

for all points (x, y) in a certain region R of the xy-plane, and for arbitrary 
values of t. In other words through each point of R there passes one and but 
one of the curves (12) in a given direction t. All six of the functions are 
supposed to be continuous and have continuous first derivatives. 

A family of curves of this sort will consist of solutions of a system of 
differential equations of the form 

(15) x, = cost, y, = sinr, t, = <£(», y, t). 

For from equation (13) and the expression for s t , the curvature t, = — is seen 

to be a function of I, u, v, and by equations (14) of x, y, r. Conversely, 
when the three equations (15) are given, it follows from the well-known 
theory of differential equations that there is but one two-parameter family 
of curves which satisfy them. The assumption that the curves (12) are the 
extremals of the integral (3) , and the assumption that the extremals must satisfy 
a given system of differential equations in the form (15), are therefore 
equivalent. 

If the extremals of the integral are to be the solutions of the system (15) , 
then from (11) the function f must satisfy the equation 

(16) /*sinT -/ J ,co8T+/ XT cosT+/ yr sinT + (f + f„)<f> = 

identically in as, y, t. This is a partial differential equation of the second 
order for f and not easily solvable by well-known methods. But if it is dif- 
ferentiated partially for t, another differential equation is found, 

(17) M x cost + My sin t + M T <j> + M<j> T = 0, 
which is linear in the function 

(18) M=f + f„ 
and its derivatives of the first order. 
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The equations (17) and (18) can now be solved by a method considerably 
simpler than that of Darboux since it avoids the use of the theory of partial 
differential equations. By means of equations (12) and (13), M(x, y, r) 
may be thought of as a function of t, u, v. From equation (17), then, 

since r s = <j> and M t = M,^xi + y\. Therefore 

(19) M = G{u t v)e- s ^^ TTf,u , 

where x, y, r in <f> 7 are supposed to be replaced by their values in terms of 
t, v, v, and the integral in the exponent stands for some particular function of 
t, u, v whose derivative for t is <f> r ^x\ + y\. The function G(u, v) is an arbi- 
trary function of u and v. The determination of M as a function of t, u, v 
in this manner means that it is also determined as a function of *, y, r on 
account of equations (14). 

It follows of course from the foregoing discussion that the function 
M{x, y, t) denned by equation (19) satisfies the differential equation (17), 
but this can also be readily verified by direct substitution. The first member 
of (17) becomes simply 

M u (u x cos t + u y sin t + u T ^) 

+ M v ( v x cos t + Vy sin t + v T <f>) 

+ M t {t x cost + t y sin t + t r <l>) + M<f> r , 

which vanishes since u and v are constant along an extremal and 

dt 
M t (t x coa t + t y sin t + t T <p) = M t -=- = — M<f> r , 

from equation (19). 

It should perhaps be mentioned that the functions u and v above are not 
the only ones which might be used in the determination of M. If « and v 
are any two independent solutions of the abbreviated equation 

(20) 4>a.COST + <I> y 8inT + & T <f> = 0, 

M 
and if ^r is a particular solution of (17), then it is true that -r- , u, v, thought 
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of as functions of the four variables x, y, t, M, are independent solutions of 
the auxiliary linear differential equation 

<t>* cos t + <S>j, sin t + & r <f> — & M M4> T = 0. 

From the theory of linear equations* it follows that the most general solution 
M(x, y, t) of equation (17) is found by solving the equation 

M 

(21) ^ = G(u, v) 

for M. The functions « and v from equations (12) , and the exponential factor 
in equation (19) are particular solutions u, v, ^ of equations (20) and (17), 
respectively. 

It remains to derive the value of / from that of M by means of the 
equation (18). This is a differential equation of the second order with re- 
spect to t, having as its most general solution t 

(22) /= /sin(T - \)M(x, y, \)dX + A(x, y) cost + B(x, y) sinr, 

where A and B are arbitrary functions of x and y. It is not certain that the 
value of /just given satisfies the Euler equation (16). In feet if the expres- 
sion (22) is substituted in (16), the equation 

(23) |^ - |? = f\M x cos X + M v sin \)dk + £(*, y, t) M (x, y, t) 

is found. The function on the right is independent of t since its derivative 
for t is the expression (17), and consequently equation (23) may be written 

< 24 > |£ - 1£ -♦<«.*<>) Jr<«.y.o). 

If Aq and B are any two particular solutions of this last equation , then it can 
readily be shown that the differences A — A , B — B„, in which (A, B) is 
any other solution, satisfy the equation 

B(A - AJ m KB - B ) 

dy dx ' 

* See for example Jordan, Covin d' Analyse, vol. 8 (1896), p. 814. 

f For a discussion of linear differential equations with constant coefficients, see Murray, 
Differential Equation*, chap. VI. 



81I1T. 
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and consequently that the most general solution of (24) is 
. . dd _ _ dd 

A=A * + Tx' B = B < + Sy' 

where 6 is an arbitrary function of * and y . 

Hence the most general function f{x, y, r) for which the integral (3) has 
a given set of extremals (12), is given by the equation 

(25) /= jf 8in(r - \)M(x, y, \)d\ + ( A + f|) C08T +(^> +§J) 

In this equation 

M(x,y, t)= G(u t v)e 

where G (w, v) is an arbitrary function ofu and v, and the arguments t, u, v 
in the expression on the right can be replaced by their values (14) in x, y, r. 
The pair of functions, A§ (x, y) and B (as, y) , are a particular solution of 
the differential equation 

and 6 is an arbitrary function of x and y. 

3. Further conditions determining / uniquely. It is natural 
to ask at this point what further conditions on the function f would suffice to 
determine it uniquely. In order to determine the function G(u, v) suppose 
that the region R has a simply closed curve 

(B) x = g(a), y = h(a) 

as its boundary, which nowhere satisfies the Euler differential equation (11). 
It is evident that no extremal can be a closed curve in Ii, otherwise after a 
complete circuit the same values of as, y, t would correspond to different 
values of t, which is contrary to the hypotheses upon the functions (12) . It will 
be supposed, furthermore, that for any fixed values of u, v corresponding 
to a point of B and a direction through that point, the functions (12) of t 
can be continued until the extremal encounters the boundary curve B. No 
extremal E can cross the curve B at a point of tangency with B. For two 
curves tangent to each other cannot cross unless they have the same curvature, 
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and by hypothesis B can never satisfy the equation (11) which defines the 
curvature r, of E. It follows then that every extremal of the set (12) crosses 
the boundary curve B without being tangent to it, at one point at least. 

Let now the values of M (x, y, r) be assigned for any point (x, y) on 
the curve B and for all values of t. The only restriction upon the values 
assigned to M, aside from the restriction that M(a, t) as a function of a and r 
shall be continuous and have continuous first derivatives, will be that the 

M 

values of the quotient --- shall be the same at all points which the extremal 

E has in common with the curve B. Thus if the region R were a circle and 

M 

the extremal E a straight line, -r would be the same for the values of x, y, r 

on the straight line at both its intersection points with the circle. This 

M 

restriction is necessary because from equation (21), — is to be everywhere a 

function of u and v alone. 

It follows now easily that for every value of « and v defining an extremal 
E in the region R, the value of G(u, v) is completely determined. For the 
extremal E intersects B at least at one point, and there the value of M has 
been assigned, so that the value of G is completely determined by equation 
(21). 

The analytical expression for Q(u, v) can be determined in an actual 
case by solving the equations 

(26) « = u[gr(o), A(a),T], v = v[g(a), A(o), t] 

for a and r and substituting the results so found in the expression 

M(a, t) 
*0(a),A(.),T]- 

The values of u, v defined by equations (26) at a point b where an ex- 
tremal E crosses B, are those belonging to E. The functional determinant 
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is different from zero at the intersection point. For along the extremal E 

u x cos t + u y sin t + u T <f> = 0, 
v x cos t + VySin t + v T <l> = 0, 
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and therefore 

cos t : sin r : <f> = 

so that the functional determinant above is 



(— cost^„ + 8inTA a )- 



This cannot vanish at a point b where E crosses B, since at such a point B 
and E are never tangent. Consequently from the theory of implicit 
functions* the solutions of the equations (26) are continuous functions of u 
and v with continuous first derivatives in the neighborhood of the values «, v 
which define any extremal E in the region JR. 

The only arbitrary part of/ besides G{u, v) is the function d(x, y), 
and the values of 6 can be determined by prescribing the transversals to any 
one-parameter family of curves C which simply cover the region Ji, that is, 
which cover the region in such a way that through each point there passes one 
and but one curve. In such a field of curves the direction angle t is a function 
of x and y, and the transversal curves in the field are by definition the curves 
of a one-parameter family whose direction angles r' satisfy the equation 



(27) 



(/cos t — / sin t) cos r' + (/sin r + f T cos r)sin t' = O.t 



If the transversals are prescribed as a one-parameter family of curves, 
w(x, y) = constant, simply covering the field and nowhere tangent to the 
original family, then ¥ is given as a function of x and y by the equations 



10, 



COST' = 



Vwi + w£ ' 



sin r = 



w u 



y/vt + »* ' 



and is different at every point from r(x, y). By substituting the value (25) 
in equation (27) it follows that the function 6 must satisfy the differential 
equation 

dff dff 

(28) j- cos t' + a- sin r'. + Q(x, y) = 0, 



dx 



*y 



* Goursat-Hedrick, A Course in Analysis, vol. 1, §26. 

f Bliss, A new form of the simplest problem of the calculus of variations, Transactions 
Amer. Math. Soc, vol. 8 (1907), p. 413. 
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where 

Q(x, y) = A^coar 1 + B sin-r' + / sin(T / — \)M(x, y, \)d\. 

The most general solution of equation (28) can readily be determined if 
a particular solution x( x > V) * s known. For the function w(x, y) is a solution 
of the abbreviated equation 

(29) _ cost' + ^ sinr' = 0, 

and w and 6 — x are two independent integrals of the auxiliary equation 

4> a .co8T / + 4> y sinT' — Q&f = 0. 
The most general solution of (28) is therefore determined by the equation 

where H is an arbitrary function of w. 

Along any particular curve C belonging to the original field, w varies 
monotonically, for on account of equation (29) it follows that 



dw 
~dl 



-j- =:w x cos t + w y sin t * 0. 



If this inequality were not true then at some point the ratios cos r : sin r and 
cos r 1 : sin ¥ would be equal and the transversal at that point would be tangent 
to C - Consequently s is also a monotonic function of w. If now the 
values of f as a function of * are given along C , then the arbitrary function 
H is determined except for an additive constant as a function of * by 
the equation 

f(s) = / sin (t - X) M(x, y,\)d\+(A 1t + x^cos t + (B + x,)«in t + -^ , 

and therefore, from what has just been said, as a function of w. 

The results of this section are as follows : 

If the, boundary of the region R consists of a simply closed curve B which 
(it no point satisfies the JSuler equation, then the arbitrary function G{u, v) 
is completely determined when the values which M —f + f„ takes at points 
of the boundary B for all values of t, are assigned. Furthermore, the func- 
tion 6{x, y) is determined, except for an additive constant, when, for any one- 
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parameter family of curves C which simply cover the region It, the family of 
transversal curves is assigned and the values off along a particular curve Cq 
are given. The function f is then uniquely determined in all that portion of 
the region R which is swept out by the transversal curves which intersect C<j. 

4. The integrals having Btraight lines as extremals. As an 

application of the results of the preceding section consider the integrals which 
have as extremals the straight lines of the plane 

x = t cos v + u sin v, y — t sin v — u cos v. 

Here u is the distance of the line from the origin, and v is the angle which it 
makes with the positive x-axis. The solutions of these equations for t, u, v 
in terms of x, y, r are 

t = x cos t + y sin t, m = x sin t — y cos t, v == t, 

and the differential equations (15) are simply 

dx dy . dr 

-j- = cos t, -£ = sin t, -=- = 0. 

as as as 

The value of M(x, y, t) from equation (19) is 

M = G(x sin t — y cos t, t)), 

the value ijr = 1 being an admissable value of yfr. Finally Aq and 2? can be 
taken equal to zero since the function 4>(x, y, 0) in (24) vanishes. The most 
general function f whose integral (3) has straight lines as extremals is 
therefore 

(30) /= / sin (t — X) G(x sin X — y cos X, \)d\ + X cos t + y sin t, 

where G is an arbitrary function of X and u = x sin X — y cos X, and is an 
arbitrary function of x and y. 

If the values of M(x, y, t) at points of a circle 

x* + y* = a* 

are assumed to be always equal to unity, then from the results of the last 
section G(u, v) must be identically equal to unity for all values of u and v 
corresponding to straight lines which cut the circle. If the straight lines 
parallel to the x-axis are assumed to have the straight lines parallel to the 
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y-axis as transversals, then in equation (28) t = 0, ¥ = ir/2, and the equation 
reduces to 

d$ _ 
dy ~ 

Consequently 6 is a function of x alone. If the values of f on the line y = 
are assumed to be all equal to unity, then from (30) 

1 = 0,, 6 = x + c, 

and it follows that the function / is everywhere equal to unity. 
The integral 



/' 



<s therefore the most general integral (3) for which : 

a) straight lines are extremals, 

b) the function M(x, y, r) equals unity for values (x, y) on a circle with 
its center at the origin and for all values of r, 

c) the family of straight lines y = c has as its family of transversals the 
lines x = c, 

d) on the line y = the values of f(x, 0, 0) are unity. 

One of the problems which Darboux solved * was the determination of 
the form which the length integral on a surface has when the geodesies on the 
surface are represented in the plane of the parameters of the surface as 
straight lines. If the variables x and y are the parameters, the length integral 
will be an integral (3) in which 



(31) f{x, y, t) = s/EcosPt + 2jFcos t sin t + G sin*T, 

E, F, and G being functions of x and y only. The problem is then to de- 
termine the functions /which are simultaneously of the forms (30) and (31). 

For (31) the value of M is 

EG-F* 

and this must be a function of t and u = x sin t — y cos t. The same must 
be true of M~ i/Z , that is, 

/n ^ N Eco&r + 2^ cost suit 4- (rsin'T _. . 

(32) ( i?a - f*)» = *<"' T >- 

•Loc. ctt., J 607. 
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The functions E> F, and O can be determined from this identity if 4>(u, t) is 
any homogeneous expression of the second degree in u, cos t, sin t. For since 
u — x sin t — y cos t, 4> may be arranged in powers of cos t and sin t, 

(33) 4>(«, t) = Aco&r + 2J3cost shit + Csin'T, 

and it follows that 



(^1(7 -#*)»' " (-4C7 - 5*)*' (^C-^ a ) 4 * 

But it is also true that <E> must be such a homogeneous expression in m, cos t, 
sin r. For by differentiating the left member of (32) for t it follows that 

and therefore 

<*u* u + 3«J*««t + iKi^-uQ^ + 4>„) + (*„ T - 3«*. T + 44> T ) = 0. 

This is an identity in x, y, t and must also be an identity in r, u, u r , since a 
set of values (x, y, t^ can be determined which gives t, «, « T any arbitrarily 
assigned values. Hence by equating coefficients of powers of u T to zero, it 
can be readily seen that 

$ = «tfT(r) + llL(r) + M{t), 

where 

dK n d*L r . d*M A7Lr A 

^T = 0, ^ + ^ = 0, -^ + 4^=0, 

from which it follows that 2T is a constant while L and 3/ are any homogen- 
eous expressions of the first and second degrees, respectively, in cos t and sin t. 

Hence the only possible length integrals for surfaces whose geodesies are to 
be represented by straight lines in the xy-plane, i. e. , the plane of the para- 
meters of the surface, are integrals of the form 



V4»(t«, cost, sin t) 
/(*, y, t) _ AC _ B i 

In this expression 4> may be any homogeneous polynomial of the second order 
in the three variables u = x sin r — y cos t, cos r, sin t, and A, B, C are the 
coefficients of its expansion (33) in powers of cost and sinr. 

Princeton University, 
Princeton, N. J. 



